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ABSTRACT 

We calculate the very high energy (sub-GeV to TeV) inverse Compton emission of 
GRB afterglows. We argue that this emission provides a powerful test of the currently 
accepted afterglow model. We focus on two processes: synchrotron self-Compton (SSC) 
emission within the afterglow blast wave, and external inverse Compton (EIC) emission 
which occurs when flare photons (produced by an internal process) pass through the 
blast wave. We show that if our current interpretations of the Swift XRT data are 
correct, there should be a canonical high energy afterglow emission light curve. Our 
predictions can be tested with high energy observatories such as GLAST, Whipple, 
H.E.S.S. and MAGIC. Under favorable conditions we expect afterglow detections in 
all these detectors. 

Key words: Gamma Rays: bursts— ISM: jets and outflows-radiation mechanisms: 
nonthermal 



1 INTRODUCTION 

EGRET detected more than 30 Gamma-Ray Bursts (GRBs) with GeV emission (Schneid et al. 1992; Sommer et al. 1994; 
Hurley et al. 1994; Schaefer et al. 1998; Gonzalez et al. 2003). The highest energy photon detected was the 18 GeV photon which 
arrived 4500 seconds after the trigger of GRB 940217 (Hurley et al. 1994). These observations motivated many interesting ideas. 
Some focused on prompt high energy photon emission, e.g., synchrotron-self-Compton (SSC) emission or inverse Compton 
scattering of photons emitted by one shell by electrons in another shell ( [Takagi fc Kobayashi 2005] ). Others focused on high 
energy afterglow processes: the interaction of ultra- relativistic protons with a dense cloud (Katz 1994), SSC in early forward 
and reverse shocks (Meszaros & Rees 1994), electromagnetic cascade of TeV 7— rays in the infrared/microwave background 
(Plaga 1995), synchrotron radiation of ultra-high energy forward shock protons (jTotani 1998[l . and inverse Compton scattering 
of prompt 7— rays by reverse shock electrons (Beloborodov 2005). 

Two kinds of high energy afterglow emission models have been discussed extensively. The first is SSC emission. Motivated 
by the successful detection of an optical flash in GRB 990123 (Akerlof et al. 1999; Sari & Piran 1999; Meszaros & Rees 1999), 
Wang, Dai & Lu (2001a, b), Pe'er & Waxman (2005) and Kobayashi et al. (2007) calculated SSC emission from the reverse 
shock. Granot & Guetta (2003) and Pe'er & Waxman (2004) applied these ideas to GRB 941017. The high energy SSC 
component of the forward component was calculated by Dermer, Chiang & Mitman (2000), Sari & Esin (2001), and Zhang 
& Meszaros (2001b). The second family of models involves the external inverse Compton (EIC) process. These include 
Comptonization of the prompt photons by the forward shock electrons (Fan, Zhang & Wei 2005b), and upscattering of far- 
UV/X-ray flare photons (assuming that they originate in internal shocks) by the forward shock (Wang, Li & Meszaros 2006; 
Fan & Piran 2006b). 

Most of the above calculations were based on the standard afterglow model. However, recently. Swift has detected nu- 
merous GRBs whose early (first 10'' s) afterglow emission cannot be reproduced within the standard model ([Meszaros 20061 
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Figure 1. Schematic cartoon of the X-ray light curve of a GRB and its afterglow, based on Swift XRT data (see Zhang et al. 2006 and 
Nousek et al. 2006 for similar plots). Also shown is a schematic optical light curve, which often does not show the same breaks as the 
X-ray light curve (Fan & Piran 2006a; Panaitescu et al. 2006; Huang et al. 2007). 

IPiran fc Fan 20071 [Zhang 2007| ) . Various modifications of the standard model have been put forward to explain the observa- 
tions. However, none are compelling and the validity of the whole model is now in question. 

High energy emission provides a new window into afterglow physics and can provide an independent test of models. 
Motivated by this, we calculate the predicted high energy afterglow emission in different scenarios. We show that there is a 
canonical high energy GRB afterglow hght curve which ought to be observed (see Fig. I15|) . The detection of the predicted 
high energy emission features by observations with GLAST or ground-based gamma-ray detectors would enable us to test the 
validity of the overall model as well as the specific modifications that have been put forward to explain Swift observations. 

The paper is structured as follows. In section 2 we review Swift GRB afterglow observations and their interpretation. 
In section |3] we describe the methods we employ for careful calculations of the inverse Compton effect; this section may be 
skipped if one is interested only in the results. In section 4 we calculate the SSC emission of the forward shock, and in section 
5 we calculate the possible high energy emission associated with X-ray fiares, including both SSC emission from within the 
flare and EIC emission from the forward shock. In section [6] we discuss the prospects for detecting high energy afterglows by 
GLAST and ground-based telescopes. We conclude in section 7 with a summary. 

2 SWIFT GRB AFTERGLOW OBSERVATIONS 

In the pre-Swift era, most of the afterglow data was collected hours after the GRB. These data were found to be consistent 
with the external forward shock model, though sometimes energy injection, a wind profile, or structured or patchy jets had to 
be invoked to account for the observations (Piran 2004). The Swift satellite has changed the situation. The X-ray Telescope 
(XRT) and the UV/optical Telescope (UVOT) onboard this satellite can slew to the direction of a GRB in real time and 
record the early broad band afterglow light curves. A schematic X-ray afterglow light curve based on the XRT data has been 
summarized by Zhang et al. (2006) and Nousek et al. (2006) (see Fig. [1} and consists of the following features: A very early 
sharp decline (phase-I); A shallow decline of the X-ray afterglow (phase-II); A "normal" decay phase (phase-Ill), possibly 
followed by a jet break (phase-IV); Energetic X-ray flares (phase-V), which may show up during any phase. Note that not 
all of these features have been detected in every burst. We focus here on the most remarkable of the new features: the slow 
decline (phase-II) and the flares (phase-V). Both are expected to have associated signatures in the high energy emission. 

In about half of the Swift GRBs, the X-ray lightcurves show an extended flattening (phase-II). In most cases, but not 
all, there is no change in the spectral slope when the light curve makes a transition from the shallow phase-II segment to 
the "normal" phase-Ill segment. The usual interpretation of the shallow phase is that it involves energy injection into the 
blast wave ( |Zhang et al. 2006, INousek et al. 20061 iGranot fc Kumar 20 06). An alternative possibility is that the parameter 
ee, which measures the fraction of shock energy transferred to the downstream electrons, varies with tim^, as would be the 

^ Note, however, that for some GRBs the break in the X-ray light curve is not accompanied by a break in the optical light curve (see 
Figure [TJ. The interpretation of this chromatic behavior is less clear. In the present work, we focus on the cases in which the X-ray and 
optical light curves break achromatically. 
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case if this parameter is shock-strength dependent. In either case the corresponding SSC emission of the forward shock would 
be different from the one anticipated in the standard afterglow model, as well as from each other. 

Energetic X-ray flares (phase-V) have been detected in several pie-Swift GRBs and in about half the Swift GRBs 
(jPiro et al. 20051 Burrows et al. 20051 IGaUi fc Piro 2006] IChincarini et al. 2007)l . The rapid decline of the flares suggests that 
they arise due to "late internal shocks" resulting from reactivation of the central engine (|Fan fc Wei 2005] |Zhang et al. 2006] 
IFan et al. 2"0 i35a King ct al. 2005", Nousek et al. 2006 Falcone et al. 2006IIPerna et al. 20 06 Proga fc Zhang 2006| |Zou, Dai fc Xu 2006[ 
IWu et al. 2007.,Dai et al. 2006..Gao fc Fan 2006..LazzatI fc Perna 2607]|Krimm et al. 2007. .Chincarini et al. 2007)l . An alter- 
native interpretation is that the X-ray flares arise due to refreshed shocks (jPiro et al. 20051 IGalh fc Piro 20061 IWuet al. 20071 
IGuetta et al. 2007]) . Once again, the GeV emission can serve to distinguish between the models. 

If the flares are produced by internal shocks, most of the up-scattered photons would arrive after the far-UV/X-ray flare. 
The high energy photons in this scenario will be produced by scattering of the flare photons in the external shock. In the 
EIC process, the duration of the high energy emission is stretched by the spherical curvature of the blast wave (Beloborodov 
2005; Wang et al. 2006; Fan & Piran 2006b) and is further extended by the highly anisotropic distribution of the up-scattered 
photons (Fan & Piran 2006b; see our Fig ll2l for a comparison). For the latter effect, see Aharonian & Atoyan (1981), Ghisellini 
et al. (1991) and Brunetti (2001) for details. 



3 SELF CONSISTENT COMPUTATION OF INVERSE COMPTON SCATTERING WITH 
KLEIN-NISHINA SUPPRESSION 

The relativistic electrons that are present in any synchrotron source will also produce very high energy photons via inverse 
Compton scattering (either SSC or EIC). We turn now to a calculation of this emission. When the energy of the electrons 
and the seed photons is sufficiently large it is necessary to take into account the Klein-Nishina correction to the scattering 
cross section. We also need to include the inverse Compton cooling in calculating the energy distribution of the relativistic 
electrons. 

The essential problem is to calculate carefully the Compton parameter, Y, the ratio between the power loss through 
inverse Compton scattering and synchroton radiation (-Fic(7e) and Ps(7e) respectively): 



y(7e) = Pi'c(7e)/Ps'(7e). 



(1) 



Throughout this work, the superscript ' indicates that the quantity is measured in the rest frame of the emitting region. In 
the regime of Thomson scattering, y is a constant, independent of the electron Lorentz factor 7e, and one obtains a constant 
reduction in the amplitude of the synchrotron emission compared to the case with no inverse Compton scattering. This makes 
computations relatively easy. However, in the general case, since the Klein-Nishina correction to Y depends on 7e, the effect 
of inverse Compton scattering on the spectrum and on the electron energy distribution is non-trivial. 
The power emitted in synchrotron radiation by an electron with Lorentz factor 7e is: 



P^M = (7e - l)<TTB"c/(67r), 

where B' is the strength of the magnetic field. The corresponding spectral energy distribution of the radiation is 



P,'(l/',7e)dl^' = P,'(7e)F 



^'s(7e 



where 1^5(7, 
F(x) = X 



= 37eeP7(47rmec), 
if5/3(C)dC, 

and K^/'j,{Q is the modified Bessel function. 

The power emitted via inverse Compton scattering is given by: 

Pi'c(7.) = 



nv. — — -—dv. 



- dtdv' 



(2) 



(3) 



(4) 



(5) 



where is the frequency of the photon after scattering. The quantity dN'^ / dtdv[^ is the scattered photon spectrum per electron 
(Blumenthal & Gould 1970). It is related to the spectral energy distribution of the inverse Compton radiation emitted by an 
electron: 



PL{vL,le)dv[c = hv[^ 



dNi, 
dtdv' 



dv[c 



(6) 



We define the auxiliary quantities g = ^^hv' / {meC?), f = hv[^/ {'yem^f?) and q = //[4(7(1 — /)], where hv is the photon 
energy before scattering. The factor g determines the regime of scattering, with the Thomson limit corresponding to g <^ 1. 
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The factor / satisfies hv' /{-f^m^c^) ^ / ^ 45/(1 + Ag) (Jones 1968; Blumenthal & Gould 1970). We can express dN'^ / dtdv[^ 
in terms of these quantities and in terms of the frequency distribution of the seed photons n^/ : 

To complete the calculation we need to know the frequency distribution of seed photons n^i . For EIC this is simple since 
the photons originate from an external source. For SSC, however, the situation is more complicated. This is because the 
photons are produced via synchrotron emission by the same electrons that are participating in inverse Compton scattering. 
The additional cooling of these electrons by IC influences their energy distribution and thus their synchrotron emission. We 
have solved this problem by two different approaches. First, we have used a simple "instantaneous" approach which involves a 
single integral equation. This method, which we describe in section 3.1, is conceptually simple and computationally fast. It is, 
however, approximate. We then describe in section 3.2 a more detailed and general dynamical approach. This more accurate 
method is the one we have used for all the calculations presented later in this paper. However, the two methods give very 
similar results in a very wide energy range, as seen in Fig|2] 



3.1 Instantaneous approximation 

In this approach we assume a functional form for the electron energy distribution n(7e) produced through acceleration in the 
shock front, and consider its instantaneous modification due to coofing. An electron of Lorentz factor 7e has a cooUng time 
given by 

If tc(7e) is longer than the dynamical time t'^ ~ R/Vc, where R is the radius of the shock front relative to the central engine 
and r is the bulk Lorentz factor of the outflow, then the electron produces synchrotron and IC emission for the entire time 
t'^. However, when tc(7e) is shorter than t'^, the electron radiates only for a time f^(7e). Thus, the total spectral radiation 
density produced by all the electrons in the fluid is given by 

/oo 

The spectral power distributions Ps'(//',7e) and Ps(u' ,-^e) are calculated as described earlier. For the inverse Compton power, 
we write eq. (|6]) as Pl^{i'[^,-ye)dv[^ ~ (1 + g)cU,jia(p' ,-ye)du[^, where a{u' ,-ye) is the Klein-Nishina cross-section, which is equal 
to 



, 3 , (l + g) 



Equation ((9l is an integral equation, since the function P(^{v','ye) inside the integral itself depends on C/^/. The quantity 
7e,min is the Smallest 7e down to which electrons are present. In dealing with equation ((9| we need to consider two cases (see 
Sari, Piran & Narayan 1998 for details and for the definitions of quantities): 

Slow Cooling: In this case, electrons with 7^ — 7^ have a cooling time tc(7,n) > t'd- Then, 7e,min = 7m, and we may use 
equation (|9} directly with 7e,min = 7m and n(7e) given by the original energy distribution produced in the shock. 

Fast Cooling: Here, all electrons with 7e ^ 7m have tc(7e) < t'^- Therefore, electrons will continue to cool below 7„i to a 
minimum 7e,miii such that 

tc(7<!,min) = t'a. (10) 

Now, for the range 7e,min ^ 7e < 7m, all the electrons are available for radiating. Initially, most of the electrons are at 7m, 
and as these electrons cool each electron will pass every 7e between 7m and 7e,min (where all these electrons accumulate). 
Hence we have 

n(7e)~n(7m), 7e,miii ^ 7e < 7m. (11) 

As usual, we assume a power-law distribution for the electron Lorentz factor: 

n(7e)d7e OC Je'^d-Je, 7e > 7m, (12) 

for which 7™ is given by (Sari et al. 1998) 
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Equation ([9]) may be solved numerically via an iterative method. The algorithm proceeds as foUows. We begin with some 
reasonable initial approximation for U^i . Using this, we compute -Pi'c(7e), t'dje) and 7e,min- Then, we compute the spectral 
distributions -Ps'(i^',7e) and P4(t^',7e) for all 7e ^ 7e,min and obtain via equation ([9} a new approximation for U^i . We take 
this U^' , or (for smoother convergence) a suitable linear combination of the new and old C/^/ , as the current approximation 
for Ui,' and repeat the steps. The iteration usually converges fairly quickly. 

This approach can be combined with any desired model for the GRB fireball and afterglow dynamics. We have used the 
dynamics described in Sari et al. (1998), except that we multiplied the calculated fluxes by a factor of 1/4 (cf., Yost et al. 
2003). 



3.2 Dynamical approach 

In this approach we follow dynamically the electron distribution as a function of time (|Moderski et al. 2000p . The main 
uncertainty is from the approximation for the initial distribution of the newly shocked electrons as a function of time. Lacking 
a better model, we assume that the electrons are accelerated at the shock wave initially to a single power law distribution; 

Q = Kj-" for 7„ ^ 7e ^ 7a/, (14) 

where the maximal Lorentz factor is given by 'yM ~ 4 x 10^ B' (Wei & Cheng 1997). The normalization factor satisfies: 
K ~ 47r(p— l)_R^nm7m~^ , and nm is the number density of the medium. We now follow the evolution of the electron distribution 
using: 

where 

[1 + V(7e)]7e - ^, (16) 



dR 67rmoc2 /3rr ' v 'wj ,c ^ 



r is the bulk Lorentz factor of the shocked medium, and /3r = y'l — l/T^, B'^ /Sir is the magnetic energy density. As usual, 
we assume that a fraction (en) of the shock energy density is converted into energy of relativistic electrons (magnetic field). 

To complete the calculations we need the location and the Lorentz factor of the blast wave as a function of time. The 
dynamics of the blast wave is obtained by solving the differential equations presented by Huang et al. (2000). The possible 
(but poorly understood) sideways expansion of the ejecta is ignored. We then calculate the electron distribution using eq. (|15p 
and the supplemental relations. The quantity n^i needed in eq.© is calculated via (for simplicity, we consider only single 
scattering): 



7M ^/ 



"-'^-rr 'f^'T / n(7e)n^)d7e, (17) 

where n-/^ ^ {AF + 3)NjJ{'iTvR^/3) « 3riV.^^/(7r_R^) (The term 4r + 3 is introduced by the shock jump condition), T' ^ 
i?/(12rc) is the time that the synchrotron radiation photons stay within the shocked medium, and 7miii ~ 3 is the Lorentz 
factor below which the synchrotron approximation becomes invalid. 

Once we know the energy distribution of the electrons, we calculate the synchrotron and inverse Compton emission, 
including synchrotron-self-absorption, and we integrate the observed flux over the "equal-arrival time surfaces" (Rees 1966; 
Waxman 1997; Sari 1998; Granot, Piran & Sari 1999). In the current code, we did not take into account the influence of 
the synchrotron-self-absorption on the electron distribution, as that done in Pe'er & Waxman (2005). However, with typical 
GRB afterglow parameters adopted in this work, for lO^sec < t < lO^sec (at late times, the high energy emission are 
usually too low to be of our interest), it is straightforward to show that the random Lorentz factor of the electrons emitting 
at the synchrotron-self-absorption frequency (IChevalier fc Li 20001 ISari fc Esin 200T|) is < 100. The modification of the low 
energy electron's distribution through the synchrotron-self-absorption is thus unlikely to influence the high energy spectrum 
signiflcantly. 

In Fig. 2, we compare the spectral distributions calculated via the simple instantaneous approach of sec. 3.1 and the 
more detailed dynamical approach of this subsection. The two methods are clearly consistent with each other. This gives us 
confidence in the validity of both calculations. Note that the multiple inverse Compton scattering are ignored in our dynamical 
approach but are included in the Instantaneous approximation. The consistence between these two approaches suggests that 
the multiple inverse Compton scattering is not important, at least for the typical GRB afterglow parameters (see also Sari & 
Esin 2001). 

In the case of EIC, the seed photon energy distribution is not influenced by the electron energy distribution. From this 
point of view the calculations are simpler. However there is another complication, viz., for the cases of interest to us, the seed 
photons are highly anisotropic in the rest frame of the blast wave. The spectrum of radiation scattered at an angle relative 
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Figure 2. Comparison of afterglow spectra calculated witli tlie instantaneous approximation (sec. 3.1, thin red lines) and the dynamical 
approach (sec. 3.2, thick blue lines). Model parameters are listed on the plot, and the timescales corresponding to each set of spectra are 
identified next to the curves. Note that all the results shown in later plots were obtained with the full dynamical approach. 

to the direction of the photon beam penetrating through this region is ( [Aharonian fc Atoyan 1981[ ): 

dtdu^^^dil' ~ I6n^i V' ^^2(1-0 6e(l-^) + 62(1-^)2^' ^ 

where dO! = 2-k sin OscdOsc, C = ^^'eic I ileiTieC^), bo = 2(1 — cos 9sc)^Ehu' /{rrieC?), cos ~ (cos ^—/3) /(I— /3 cos 9),6\s the angle 
between the line of sight and the emitting point, (5 is the velocity of the emitting point, and hv' <C hv'^^^ ^ 7emeC^ 69/(1 + 69). 
As expected, on integration over 9^^, eq. (|18p reduces to eq. ((7}. The energy loss rate of the hot electron beam can be estimated 
by eq. ([5J and F(7e) is governed by eqs. (|H7|) for a given n„' (see i|5.2l for details). 



4 HIGH ENERGY SSC AFTERGLOW 

The dominant source of long lasting high energy GRB afterglow emission is SSC of the hot electrons in the forward external 
shock. At early stages when the cooling of most electrons is important, the luminosity of the SSC emission, I/ggc , is related 
to the luminosity of the synchrotron radiation, i^y,,, as l|Sari fc Esin 200T|| : 

Lssc-yLsyu, (19) 

where y is the Compton parameter. The X-ray luminosity Lx is a small fraction of L^^^^ but we can use it as a proxy for the 
total luminosity. To do so we define a factor ex such that Lx = ex^syn and: 

L^sc-yLx/e^. (20) 

As long as ex does not vary significantly with time, we expect the broad-band SSC afterglow light curve and the X-ray light 
curve to have a similar temporal behavior. We expect, therefore, that Lx and Lssc should be highly correlated. This is, of 
course, confirmed by more detailed analysis, as shown below in eq. (|3ip . 

The light curve depends on the dynamics of the blast wave and in particular on the evolution with time of -£/^[„, the 
power given to the shocked electrons (see eq. I22p . We consider first the evolution expected in the standard afterglow model 
and then discuss various modifications to the model. 



4.1 Analytic Considerations 

We begin with the standard afterglow. We consider a circumburst medium with a number density profile = n*i?~'°, 
^ < 3; here, = corresponds to a constant density ISM, and = 2 to a standard stellar wind (Meszaros, Rees & Wijers 
1998; Dai & Lu 1998b; Chevalier & Li 2000), though k ~ 1.5 is still possible, as found in some supernovae (jWeiler et al. 2002P 
and in GRB 991208 (|Dai fc Gou 200ip . The quantity n, is the number density at a distance R — 1: 

= 35 . 3 (21) 

1 3.0 X 10^^ A* cm-^ for k=2, ^ ' 
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where A, — [M/10 ^Mq yr ^] [uu,/(10*cm s ^)], M is the mass loss rate of the progenitor, v^n is the velocity of the stellar 
wind (|Chevalier fc Li 2000|) . 

Following the standard afterglow model, we assume that the dynamical evolution of the ejecta has a Blandford-McKee 
self-similar profile l|Blandford fc McKee 1976]| . The power given to the freshly shocked electrons, in the blast wave is: 

_ p 1 r 7.5X 10^« ergs-i, for fc = 0; 
L ^ 6.,_ii5.,53t3 I 5xl0«ergs-i, for fc = 2. ^^^> 

where -Efc is the equivalent isotropic energy of the ejecta. We note that L^j^ depends only weakly on the density profile. Here 
and throughout this text, the convention Qx ~ has been adopted in cgs units. 

The SSC luminosity can be estimated as: 

-'^ssc ~ '^high-^'cin- (23) 

All the physics in this equation is, of course, hidden in the factor e^^.^^^, which depends, in turn, on the synchrotron and cooling 
frequencies, Vm and Vc (jSari et al. 1998|) . and on the power law index of the electron distribution, p: 



high 



^vy/i^+y), (24) 



where (jSari et al. 19961 ISari fc Esin 200l|l : 

n = min{l,(z.„K)(*'-2)/2}, (25) 



y ~ i-l + ^l + 4fjee/eB)/2, (26) 

77 = min{l,(i.„/^,)(''-2'/2}, (27) 

i^c = {l+yffc. (28) 

The ratio of the synchrotron and cooling frequencies satisfies l|Sari et al. 19981 lYost et al."2003|) : 



0.0024C;el_ie|;__2n£;fc,52f3 for fc = 0, 
0.l2C^ei_^el^_2Al_it-^ for fc = 2, 



(29) 



where Cp = 13(p — 2)/[3(p — 1)]. Note that, in all analytical relations, the time and the frequency are measured in the burst's 
frame, that is we ignore cosmological (l+z) corrections. Numerical results are presented for a canonical burst at 2 = 1. 

The X-ray band is typically above max.{i'm,Vc}- In this case the forward shock X-ray emission can be related to the 
kinetic energy of the forward shock (|Kumar 20001 IFreedman fc Wajcman 200TllFan fc Piran 2006a|) : 



(p-2)/4 p-i . .-1 (p+2)/4 (2-3p)/4 I 8.8x10 ergs s , forfc^Q; 

^..-i(l + y) K.sz i3 <; X io« ergs s-\ for fc = 2. ^^^> 



We thus have 

-^ssc 2-p ^{2-p)/4p{2-p)/4 3(p-2)/4 J 2, for fc = 0; 

^'ny^.,-l^B,-2 ^fe.53 *3 <| 1^ forfc = 2, ^"^^^ 

For a universal p ~ 2.1 — 2.3, L^^^/Lx is sensitive only to r] and 3^ and it is only weakly dependent on other parameters. 
At early times, when the cooling of electrons is important, igg^ oc t^^'^^^'^ Lx (note that in the standard afterglow model, 
Ek, n, , ee, es are all constant). Therefore a wide band SSC light curve will have a temporal behavior quite similar to that 
of the X-rays. 

Roughly speaking, the energy of the SSC emission peaks at a frequency ~ max{j^^^'^, v^^''^}, where v^''^ ~ and 
j/f « where 7c is the cooling Lorentz factor of shocked electrons. Following the standard treatment (jSari et al. 19981 

IChevaher fc Li 2000|) . we have 

^ssc ^ ^o2i ^4^4 ^1/2 r 6.2 n-''^El%tf'\ for fc = 0; 



.f^^io^^Hz(i+^)-Vi2 : \ ' r: m 



.SSC ^ n ^ ^A-K-y- / ^ nr^"E-l^t;'^\ for fc = 0; 

1.5 A-!./i'£;fe,53t3, for fc 

Note that oc n"^/"* or oc So n ~ 10^ cm'^ or A, ~ a few will shift v^^'^ to the X-ray/UV/optical band and in 

this case the SSC emission will influence the X-ray observations. It may even cause a flattening of the X-ray light curve due to 
the emergence of this new component. An example of such a case is shown in Fig|3] One can see a flat X-ray segment, which 
is rather similar to that detected by Swift. However, it is not clear that this can account for the Swift observations because in 
this case the X-Ray spectrum would vary with time (see the insert of Fig[3} . Such variations are not seen in the Swift data. 
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Figure 3. The XRT light curve in a dense ISM. The thick solid line is the predicted X-ray light curve, including both the synchrotron 
and the SSC components of the forward shock. The insert shows the corresponding X-ray spectra at three different times, as marked in 
the plot. 



The (adiabatic) standard afterglow model assumes that (i) the outflow energy is a constant and (ii) the shock param- 
eters are constant. As mentioned earlier this model is inconsistent with the shallow decline phase (phase-II). One possi- 
bility is that one of these two basic assumptions should be revised ( [Zhang et al. 2006] IFan fc Piran 2006al lloka et al. 20061 
IPanaitescu et alT 2006; No usek et al. 20061 IGranot et al. 2006|) . We consider energy injection of the form oc ([Cohen fc Piran 19991 
[Zhang fc Meszaros 2001ai ), where q = 1 represents no energy injection and q = corresponds to a pulsar/magnetar-like energy 
injection (Dai & Lu 1998; Zhang & Meszaros 2001a; Dai 2004; Fan & Xu 2006). Other q values are possible for an energy 
injection that results from slower material progressively catching up (Rees & Meszaros 1998; Kumar & Piran 2000; Sari & 
Meszaros 2000; Granot & Kumar 2006) or if an energy injection is caused by the fall-back of the envelope of the massive star 
(MacFadyen, Woosley & Heger 2001; Zhang, Woosley & Heger 2007). We also explore the situation where the equipartition 
parameters, te and es, are shock-strength dependent (i.e., time dependent though the underlying physics is far from clear 
([Piran fc Fan 2007)) . Instead of exploring the possible physical processes that lead to such a phenomenon, we simply take 
(ee, es) oc (f^, t''). These modifications lead to: 



It is straightforward to show that 

d , (6-fc)(l-t|) + (5fc-18) 



,ssc 



oc t'^"^^'^ "Hi-") ■ ^ (36) 

7^, 6(fc-l) (7fc-10) 

(i-k)\l+y)-*, (37) 



4-fe 2(4- 



Equation ([35p is one of our main results. As expected, with significant energy injection, or either increasing with time, or 
both, Lgg;, (general) is flattened. An eg decreasing (increasing) with time will also flatten (steepen) the high energy emission 
light curve. However, such a modification seems to be small and it cannot give rise to either the observed shallow decline 
phase of the x-ray light curve or to a detectable signature in the high energy component. Therefore, we focus on models with 
either time-dependent Ek or te. 

The shallow decline seen in the X-ray light curve during phase-II (Fig. 1) requires q ~ 0.5 or c ~ 0.4. In general, for 
Lx oc t"" (a < 1) we need (in the energy injection case) q = [4(a-f l)-2p]/(p-f 2) which yields L^^^ oc t-I4(«+i)-2pl/(p+2) ^ ^-a 
for p ~ 2. The high energy decline is quite similar to the decline of the X-rays. For a varying (with no energy injection, 
i.e., q = 1), we need c = (3p - 2 - 4a)/[4(p - 1)], which in turn results in Lssc oc t{p+2-i2a)/[8(p-i)l ^ ^{i-3q)/2 p ^ 2. 
The high energy decline is slightly slower than that of the X-rays in this case. 



2 One may speculate that the energy distribution index of the accelerated electrons p is also time-evolving. However, this is not seen in 
the data as the spectrum does not vary during this phase. 



9 




Figure 4. SSC radiation from the forward shock for the case of a constant density external ISM. (a) SSC light curves in the energy 
ranges 20MeV-300GeV and 0.2 keV -100 TeV, respectively, (b) Spectra at three selected times; thin and thick lines correspond to the 
pure synchrotron spectrum and SSC+synchrotron spectrum, respectively, and solid, dashed and dotted lines are at 2 X lO'^, 2 X 10* and 
2 X 10*= s after the burst. 



At least in principle, one could combine IR/optical/UV/X-ray and high energy observations to distinguish between the 
two modifications described above. For example, we have Um oc e^E^^^ and Vm'^ oc 

^4£(6-fe)/[2(4-fc)i^ If the early X-ray 

flattening was caused by (x t", Vrn and Vm'^ will decline much more slowly than in the energy injection case Eu oc The 
wide energy range of LAT onboard GLAST (20MeV-300GeV) might enable us to observe the variations of vf^'^ with time. 

It is interesting to note in passing that these modifications provide a possible explanation for some long-term puzzles in 
GRB 940217. The long-lasting MeV to GeV afterglow emission of GRB 940217 Purley et al. 1994| ) showed two remarkable 
features: (a) the count rate of high energy photons was almost a constant (b) the typical energy of these photons was nearly 
unchanged. These two features can be reproduced with c ~ g ~ 1/2 and d = k = Q (Wei & Fan 2007). 



4.2 Numerical Results 

We turn now to numerical computations of the high energy light curves. We consider, first, the standard afterglow model 
using typical parameters that seem to fit the average late afterglow: Ek = 10^^ erg, p = 2.3, te ~ 0.1, e_B = 0.003 and Oj = 0.1. 
We consider a typical burst at 2 = 1. Figures |4] and [5] depict the calculated light curves and spectra for two models of the 
external medium: a uniform density ISM and a stellar wind. In both figures, panel (a) shows the SSC emission afterglow light 
curve and panel (b) shows the spectrum. 

We consider now an energy injection model where the energy injection has the form: 

dEinj/dt = 5 X 10''^(t/100s)"°-'^ erg s"^ (38) 

for 10^ s < t < lO** s, which corresponds to q = 0.5 and Ek = 10^^ erg. Apart from q and Ek all other parameters are similar 
to those used in the standard case above. The total integrated energy injected is equal to 9 x lO^^erg ^ Ek- The resulting 
light curve is shown in Fig. [6l The SSC light curve is fiattened when the energy injection is strong enough to suppress the 
deceleration of the outflow. The numerical light curve has Lssc oc which is consistent with our analytic estimate 

Lssc oc t~°-^ for c = and q = 0.5 (see eq. ((35|l '). 

We turn now to a time-evolving shock parameter e^, and consider te varying as As shown in eq. (|35|l and in Fig. [71 
an increase with time of ee flattens the high energy emission light curve. The very small at early time not only lowers the 




Figure 5. SSC radiation from the forward shock for the case when the external medium corresponds to the wind from the progenitor 
star (A; = 2). The hne styles are the same as in Fig. |4] 




Figure 6. SSC radiation from the forward shock for the case when energy from the central engine is injected over a period of time: 
dEi^j/dt = 5 X 10''^(t/100s)~0'^ erg for 10^ s < t < lO"* s. Note that the SSC emission light curve flattens as a result of the energy 
injection. 



fraction of the shock energy given to the fresh electrons but it also suppresses the SSC emission. The resulting t~°'^ decline 
depicted in Fig. [7]is consistent with the analytic estimate for c = 0.4 (and q = 1). 

To check the consistency of the numerical and analytic results, we plot the two estimates of I/ggQ (using eq. I23p in Fig. 
[SI The analytic results (the thick lines) are a factor of 2 — 4 times larger than the corresponding numerical results (the thin 
lines). This is reasonable as some important corrections, such as the integration of the emission over "equal-arrival surfaces", 
have been ignored in the analytic formulae. 



11 




10' 10' 10' 10* 10' 10° 

t (sec) 



Figure 7. SSC radiation from the forward sliock for the case when the electron energy parameter te varies with time. The soUd and 
dashed Unes correspond to the emission in the energy range 20MeV-300GeV, while the dotted and dot-dot-dashed lines are for the 
emission in the energy range 0.2 keV-100 TeV. The shock parameters are e_g = 0.003, = 0.017 for t < 100 s, = 0.017(t/10)0 -* for 
t < lO'* s after which it saturates. Other parameters are Ek = 10^^ erg, z = 1, 9j = 0.1, p = 2.3. The parameters corresponding to the 
external mdium are marked on the plot. 




Figure 8. Comparison of numerical and analytical results for ISM. The thin lines are our numerical SSC light curves for the energy 
range 0.2keV — lOOTeV, while the thick lines are the corresponding analytical results. The solid and dashed lines are the results in the 
cases of energy injection and increasing with time, respectively. The parameters are the same as in Fig. |6] and in Fig. [7] respectively. 
For the wind medium, the results are rather similar. 



5 HIGH ENERGY EMISSION ASSOCIATED WITH X-RAY FLARES 

We turn now to GeV flares that might arise from inverse Compton scattering of the radiation associated with X-ray (or UV) 
flares. Aflhough X-ray flares (phase-V in Fig[T]) were detected even before Swift, their frequency became clear only after Swift 
began its observations. By now it is known that flares are quite common and can appear at all phases of the afterglow. At 
times the energy emitted in a flare can be fairly large. There are two main ideas to explain the origin of these flares: (i) "Late 
internal shocks" (Fan & Wei 2005, Burrow et al. 2005, Zhang et al. 2006) associated with a long-lived central engine, and (ii) 
"refreshed shocks" IjGuetta et al. 2007|l when late shells encounter the external shock and lead to brightening. 

Inverse Compton scattering of photons from an X-ray flare are possible via two distinct mechanisms. It could be the 
result of SSC emission from the same electrons that produce the X-ray flares. If the X-ray flare is produced by late internal 
shocks, then an additional source of high energy radiation is possible, viz., EIC scattering of flare photons by hot electrons in 
the external shock. We consider both possibilities. 
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5.1 SSC flares 

SSC within the same shock that produces the X-ray flare will give a high energy flare simultaneously with the low energy 
flare. This would arise if the X-ray flare results from either a late internal shock or from a refreshed shock within the forward 
external shock. 

It is difficult to predict the expected SSC emission as we have no robust estimate of the typical Lorentz factor of 
the shocked electrons that produce the X-ray flare. A critical factor is the location of the shock, which determines the 
various parameter within the emitting region. For prompt 7— rays from an internal shock, the typical radius of the shock 
is -Rprompt ~ 10^^ — lO^^cm l|Piran 19991 IPiran 200411 . If flares are produced by "late internal shocks", iJflare ~ 10^* cm is 
possible UFan fc Wei 2005)l . whereas with "refreshed external shocks", -Rflarc may be as large as 10^ cm UGalli fc Piro 20071 
IWu et al. 20071 IGuetta et al. 2007|) . 

Assuming that the soft X-ray flares are powered by the synchrotron radiation of the shocked electrons we can estimate 
the typical Lorentz factor of the electrons, 7c, m- The magnetic field, B, at -Rflarc can be estimated by: 

B ~ [2eLx/ir^RLrcc)]'^' ~ 250 Gauss e'^' Ll%r-' R^L^ir, (39) 

where e = es/ee. For this value of the magnetic field, the peak energy of the fiare photons will be at Ep ~ 0.2 keV if the 
typical electron Lorentz factor is 

7c,n. ~ 800 e-'^'j:x!/9<aLi5(^p/0-2 keV)^''^ (40) 
The energy of a typical inverse Compton photon is then 

/li^r ~ '2'ylnMp ~ 0.3GeV e-^^''L-]i^Rn,rc,i5{Ep/0.2 keVf. (41) 

Thus, high energy emission simultaneous with the X-ray fiare is expected if the emitting region is not significantly magnetized. 

Roughly speaking, the total fluence of the SSC emission of the flare shock is comparable to that of the X-ray emission, 
typically 10"'^ ~ 10"'^ erg cm ^ integrated over the aftgerglow. In late internal shocks (i.e., -Rflarc ~ 10^^ cm), a GeV flash 
accompanying the X-ray flare is possible (Wei, Yan & Fan 2006). This problem was also discussed by Wang et al. (2006), who 
assumed that 7c, m ~ 100 and obtained hUp"'' ~ 10 MeV, which they considered as uninteresting. However, as shown above, 
7e,m can be large to ~ 1000 and hVp^'^ is two orders of magnitude larger. In refreshed external shocks, a GeV-TeV flash is 
predicted because of its very large -Rflarc (~ 10^^ cm), as indicated in eq. (|4ip . More detailed analysis can be found in Galli 
& Piro (2007). 

A subtle issue that has to be checked is whether the high energy photons will be absorbed by pair production on the high 
energy tail of the flare. The pair production optical depth for photons with energy -Ecut (absorbed by the flare photons with 
energy Ea ~ 2(rmeC^) V-^cut ~ 0.5MeV ri.5(-Bcut/lGeV)"^) can be estimated as (e.g., Svensson 1987) 

llaTA^>-E„ 2 p-2 p „ ,A+, n2 ( Ep \ftlaro-lr-2ftlarol- ^''^^ ^49^ 

,^^^^^-4x10 i?,_^^Ffl_,_s.35tiD^,,3,34(^^^) (^^) , (42) 

where N^e^ = ^/a^'""^ ('1^)^""° '^^^ L^as^r^^st .^^ total flare photon number of one pulse satisfying hv > Ea, where St is the 
timescale of the flare pulse and the high-energy power-law index /3flarc ~ 1.2 has been used to get the numerical coefficient. 
Clearly, for -Rflarc ~ lO^'^cm, i.e., the refreshed shock case, the tens GeV high energy photon emission will not be absorbed by 
the flare photons. For -Rflarc ~ 10^^ cm, i.e., the late internal shock case, the small optical depth will not affect the sub-GeV 
flux unless 5ti > 25ri^gn''"=7?Lrc,i5- 



5.2 Extended EIC plateau 

We turn now to the scenario in which the X-ray flares are produced by late internal shocks (Fan & Wei, 2005, Burrow et 
al. 2005, Zhang et al. 2006). We calculate the inverse Compton scattering of these seed photons by hot electrons accelerated 
within the external shock. We assume that the X-ray flares are accompanied by far-UV emission and calculate the upscattering 
of these photons as well. A central ingredient of this scenario is that in the rest frame of the blast wave, the seed photons are 
highly beamed. We take care of this effect, following the analysis of Aharonian & Atoyan (1981). 

If the EIC emission is simultaneous with the X-ray flare (i.e., the duration of the EIC emission has not been extended 
signiflcantly) , the EIC luminosity can be estimated by eq. (I23|) . However, in the rest frame of the shocked material, the EIC 
emission has a maximum at 9sc = tt and it vanishes for small scattering angles ( [Aharonian fc Atoyan 1981 [ IBrunetti 2001 P . 
This effect lowers the high energy flux in two ways. First, a fraction of the total energy is emitted out of our line of sight 
and thus the received power is depressed (relative to the isotropic seed photon case). This yields a correction by a factor of 
2 (which we ignore henceforth). Second and more important, the strongest emission is fronti 61 ~ 1/r. Thus the peak time of 
the high energy EIC emission is estimated to be (Fan & Piran 2006b; Wang & Meszaros 2006): 



^ This could be more easily understood in the Thomson regime. As shown in eq. (43) of Brunetti (2001), in the local frame of the 
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Tp ~ R/{2V'c) ~ (4 - k)tu (43) 

where U is the time when the X-ray flare ceases. Tp, which is also proportional to the duration of the high energy peak, could 
be much longer than AT, the duration of the soft X-ray flare. 

The luminosity of the high energy flare would be lower than the simple estimate by the ratio of the durations: 



L 



(44) 



'=^<= (Tp/AT) ■ 

At 100-1000 s after the burst, the forward shock emission peaks in the far-UV to soft X-ray band, and the corresponding SSC 
emission peaks in sub-GeV to GeV energy range. A comparison of the SSC luminosity of the forward shock after but around 
tf, Lgsc (eq. (|23|l ). with L^^^ shows that the SSC emission would be stronger than the EIC emission and the wide ETC flare 
would be undetectable. 

However, if the forward shock electrons are in the slow cooling regime before the X-ray flare, their SSC emission is 
weak and the EIC flare might be detectable. In this case the total energy available for extraction in the EIC process ~ 
LeinAT -I- A^eF min{7c, 7m}meC^ is much larger than ~ LeinAT, where A^e is the total number of electrons swept by the 
forward shock at the time ~ t/ — AT and at the same time igsc much smaller than Lein. Though L^^^ may still outshine 
Lj-ic at t ~ t/, since it decreases rapidly with time (steeper than t~^, as both rj and y/{l +y) axe decreasing with time, see 
eq. (|23|l l. the EIC high energy emission may still dominate at late times. 

If the EIC emission dominates over the SSC emission, the high energy light curve will flatten, as we show below (e.g.. 
Fig. Such a flattening could arise also as a result of energy injection or due to an increasing e^- However, as we argued in 
the last section, in those two scenarios, the X-ray and the high energy emission light curves are quite similar and flattening 
should be apparent also in the X-ray signal. The EIC emission should, on the other hand, show an X-ray flare preceding high 
energy emission and not accompanying a flat X-ray light curve. 

As an example we consider the giant flare of GRB 050502b HBurrows et al. 20051 IFalcone et al. 2006|l and examine the 
expected external IC emission that will arise from such a flare. The flux of the flare, in the 0.2—10 keV energy band, can be ap- 
proximated as a steep rise: Fame « 5 x 10~^ erg s~^ cm~^ (f/680s)^ for 300s < t < 680s, a constant plateau lasting until ~ 800s 
and a subsequent sharp decline which might be due to a curvature emission component ( [Fenimore, Madras fc Nayakshin 1996| 
IKumar fc Panaitescu 2000l [Liang et al. 2 006). To calculate the EIC emission we need (see eq. [7| n'^, the distribution of the 
seed photons in the rest frame of the shocked medium. 

If the flare originates from activity of the central engine l|Fan fc Wei 2005] [Zhang et al. 200"6| one might expect that the 
radiation process is similar to that of the prompt emission. Lacking exact information on the spectrum of the flare and in 
particular on its peak energy we assume that it has a typical Band function ijBand et al. 1993)) : 

„ ^^f (T^)-<^+"-"'exp(-^'K), for^'^B^;; 

I (|^)''"-"-""""exp(-/3fl,,e-afl,,e)(^)-(i+A--), forBz.;^;.'; 

where the high-energy power-law index /3flarc ~ const ~ 1.2 and the low-energy power-law index asare ~ constant ~ 0. As the 
peak energy is not known we consider three representative values: Ep — 0.02, 0.2, 2keV. 
For a given i/'p, the parameter A in eq. (|45p is obtained from the observed flux: 

lookcv/hr „2 p 

iljl ^L-f^ flare / Ar'\ 

n,,hv dv ^ pa 2 , (46) 

0.02keV/hr ^ 

where we used lower and upper limits on hv' of 20eV/r and lOOkeV/F. This is because the self-absorption frequency is likely 
to be in the UV band (|Fan fc Wei 2005|l and the emission in the hard X-ray band is unknown for nearly all flares. 

The redshift of GRB 050502b is unknown. We assume the canonical value ai z — \ for which i^k ~ lO^'^ erg because the 
7— ray fluence is ~ lO"'' erg cm~'^ ([Burrows et al. 2005|l . For the other parameters we take n — \ cm~'^, p = 2.3, ee — 0.1, 
es = 0.01, and Qj = 0.1. Figures l9l and [TOl depict the electron distributions and the Compton parameters as functions of 7e, 
for Ep ~ 0.2 keV. The cooling effect of the X-ray flare photons on the blast wave electrons is seen clearly in these flgures. One 
sees that the energy of the electrons is depressed between 400 and 700 s and then it increases at 900 s when the cooling effect 
due to the flare photons ceases. As expected the higher the flare luminosity, the stronger the EIC cooling. Electrons with 
7e < 10® lose most of their energy via the EIC process (see the large values of the Compton parameter for these electrons in 
Figs.[9l;b) andllOjb)). 

The resulting high energy emission is shown in Fig. 1111 The SSC emission decreases during and after the flare as the 
electrons are cooled by the EIC process. Also at a later time we get contributions to the observed spectrum from higher 

shocked medium, the emissivity is proportional to (1 — cos Ssc)'^^'''''^, where (5 = p or p + 1, depending on the cooling of the electrons. 
The observed emission from an angle d is thus oc [r(l - /3 cos 6*)] sin 5(1 - cos 650)*^+'^'''^ oc -|- r2e2)-(7+'5)/2 = jF since 

cos6»sc = (cos6» - - /3cos6»). The requirement that dTldB = yields 9 « ^(2-|-(5)/5/r ~ l/F for 5 ~ 2 - 3. 
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Figure 9. The electron distribution A^^^, and the Compton parameter V(7e) as functions of the electron Lorentz factor 7c. The flare 
photons are assumed to scatter off electrons in the forward shock for the case of a uniform external ISM. Different lines represent different 
times after the initial burst, as indicated in the figure; the times are determined by dt = (1 + z)dR/ {2T^ c). The parameters of the model 
are: E]^ = 10^^ erg, n = 1 cm"'', z = 1, p = 2.3, te = 0.1 and = 0.01. The parameters of the flare are described in the text and the 
peak energy of the flare emission is taken to be Ep ^ 0.2 keV. 




Figure 10. The electron distribution A'^^ and the Compton parameter Y(7e) as functions of the electron Lorentz factor 7c. Here the 
flare photons are assumed to scatter off electrons in the forward shock for the case when the external medium is due to a stellar wind 
with A-t =0.1. Other parameters are the same as in Fig. |9] 
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Figure 11. High energy light curve arising from flare photons scattering ofi' the forward shock for the ISM/wind case (upper/lower 
panel), respectively. The parameters are the same as in Figs.[9land llOI except for Ep which takes the values marked in the figure. 



latitude regions from which the emission is weaker. The EIC emission is not simultaneous with the X-ray flare. It peaks at 
~ {4: — k)tf (see eq. (HSJ), and it lasts much longer than the X-ray flare. This temporal behavior is determined by the geometry 
of the emitting surface, the radiation spectrum and the highly anisotropic EIC emission. The lagging behavior is unique and 
if it is observed, i.e., if it is not hidden by SSC emission, it would demonstrate that X-ray flares are produced by internal 
shocks. Note that without the anisotropic correction, the EIC light curve is higher and narrower and the peak EIC emission 
is overestimated by one order of magnitude (see Fig[T2}. 

For most soft X-ray flares, the peak emission energy seems to be below 0.2 keV, i.e., they may be intrinsically far-UV 
flares. The up-scattering of the far-UV photons in the external blast wave results in strong sub-GeV emission. Fig. [13] depicts 
the resulting EIC spectrum (time integral) for different values of Ep = (0.02, 0.2, 2) keV — other parameters, including the 
luminosity of the flare in the 0.2-lOkeV band are taken to be the same. For a far-UV flare {Ep ^ 0.2 keV) the seed photons 
are much more numerous than for an X-ray flare. Consequently the resulting sub-GeV photons are much more numerous than 
those resulting from a keV flare. Therefore the EIC emission following a UV flare will be easier to detect. 



6 DETECTABILITY OF HIGH ENERGY EMISSION IN THE AFTERGLOW 

We turn now to the key question: Are the GeV to TeV high energy signals predicted by our models observable with current 
or soon to be commissioned detectors? 

Using the calculated high energy spectrum F^{t) as a function of time for any given model, we can estimate the total 
number, A'^dot, of detectable high energy photons, 

A^dct = / / ^j^Sd^t{y)dtdu, (47) 

where ti,E are the times when the observations begin and end, respectively, hvd — huu is the energy range of the detector, 
and Sdet{v) is the effective area of the detector as a function of For LAT onboard GLAST, we approximate S'dct(j^) as (see 



http:/ /www-glast. slac.stanford.edu /software/IS/glast_ lat_performance.htm): 




Figure 12. The EIC emission with (the dotted line) and without (the thin dashed line) anisotropic correction. Other lines and the 
parameters are the same as in the upper panel of Figs. [TTI except for Ep that marked in the figure. 





Figure 13. The EIC spectrum resulting from upscattering of flare photons by forward shock electrons. The parameters used are the 
same as in Fig. [9] and Fig. 1101 except for Ep which takes the values marked in the figure. 



J 500 cm^ (/i!//20MeV), for hu < 400MeV; 
^-dcU^^) - I ^g4^^2^ f^^. ^ 400MeV. ^^^^ 

We consider first the high energy SSC emission in the afterglow, which we estimated in sec. 4. For the models presented 
in Fig. |4]- Fig. [T] we use ti ~ 100 s; at earlier time the high energy emission may be dominated by the synchrotron and/or 




Figure 14. The integral of vFi, in the time interval 100 — 4 X 10* s after the burst. The solid, dotted and dashed lines correspond to 
the standard afterglow model (Figs. |4] and [Sj, the energy injection model (Fig. |6]l and the increasing te model (Fig. [7j, respectively. 



Table 1. Expected signal for GLAST from SSC emission of a GRB forward shock, where the absorption of the very high energy photons 
by IR background is ignored. These values we calculate for a typical burst with E]^ ~ 10^"^ erg (at the end of the X-ray shallow decline) 
and 2 = 1, correspondingly to the burst with a 7— ray fiuence of ~ 10~^ erg cm~^. 
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SSC emission of the internal shocks ( |Gupta fc Zhang 2007| . We choose an upper hmit of ^ 4 x 10'' s; after this time the 
SSC emission is usuaUy too low to be of interest. 

Figure [14] shows the integrated flux expected for the various SSC scenarios discussed in sec. 4, and Table [T] summarizes 
the expected number of photons that would be detected by LAT from a burst with standard parameters (see Figs. |4]- [7]) at 
2 = 1. Typically, one expects to detect a few photons above 20 MeV and very few high energy photons above 100 GeV. 

Not surprisingly, the modifled afterglow models that account for the shallow X-ray light curve in phase-II give fewer 
counts than the standard afterglow model. The reduced X-ray flux in these models (needed to explain the shallow light curve) 
causes a corresponding reduction in the high energy flux. However, we still expect a weak detection by GLAST. Such weak 
signals, of course, cannot play an important role on distinguishing between the different models. But for some extremely 
bright events, e.g., GRB 940217, the high energy observation may pose a tight constraint on the underlying physical process 
(e.g., Wei & Fan 2007). 

Considering next the high energy emission associated with flares, the time-integrated uF^ of the high energy EIC compo- 
nent is shown in Fig. [13] In the case of a uniform external ISM, for Ej, = (0.02, 0.2, 2) keV, A''dot(> 20MeV) = (0.6, 0.5, 0.2) 
and Afdot(> lOOGeV) = (3.2, 3.2, 2.6) x lO"'*, respectively. The EIC high energy afterglow component is more easily de- 
tected if the flare has a signiflcant UV component. Note that Fan & Piran (2006b) used a larger effective detection area 
(Sdct ~ SOOOcm^ in the energy range of 20MeV — SOOGeV). This overestimates Sdot for hv < lOOMeV where most of the up- 
scattered photons are expected (see eq. [48]). In the case of an external shock in a stellar wind we have for Ep = (0.02, 0.2, 2) 
keV, iVdot(> 20MeV) = (0.3, 0.3, 0.2) and iVdct(> lOOGeV) = (1.7, 1.5, 1.0) x 10"", respectively. Now the scattered far-UV 
photons are in the sub-MeV band (Fan & Piran 2006b). Therefore, the far-UV component does not increase the detected 
signal. 

As long as the flare outflow is just weakly or even not magnetized and -Rprompt 5S 10^* cm, the GeV SSC emission fluence is 
expected to be comparable to the fluence of the keV flare, typically 10~^ — 10~^ erg cm~^. With such a fluence, the GeV flashes 
(SSC emission) accompanying bright flares may be detectable by GLAST. If the flare is produced by a late internal shock we 
expect that the typical SSC photon energy is about 300 MeV. At this energy a fluence ~ 2.4 x lO'^/Sdot ~ 3 x lO"'^ erg cm 
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Figure 15. Summary of the results. The expected high energy afterglow signatures are shown in the lower panel, corresponding to the 
schematic X-ray afterglow light curve shown in the upper panel. Note that the EIC emission light curve could be outshined by the SSC 
emission of the forward shock. The SSC emission of the X-ray flares might be weak if the emission region is significantly magnetized or 
Kflare is much Smaller than 10^* cm. 

corresponds to a detection of 5 photons. So the SSC emission of a very bright X-ray fiare with a fluence ~ 10~^ erg cm~^ 
should be detected, provided that the Compton parameter is unity or larger. If the flare is produced by a refreshed shock the 
typical photon energy would be higher, up to tens of GeV. The number of these high energy photons would then be much 
smaller than in the case of late internal shocks. As a result, it might not be detectable by GLAST. 

So far we have focused on the detectability of high energy emission by GLAST. However, there are also other detectors. 
MAGIcfl, WhipplfHand H.E.S.S.EI are high energy telescopes operating at energies above 100 GeV. These Cerenkov detectors 
have very large effective areas ~ 10* — lO^m^. The expected fluxes of very high energy (> lOOGeV) photons from bursts at 
2 « 1 should correspond to the detection of 10 — 100 photons. However, this estimate ignores the absorption of the high 
energy photons by the IR-background ( Nikoshov 1962p . Given that the optical depth for a 100 GeV photon from z = 3 is > 5 
l|Primack et al. 200511 . we expect that for most Swift hursts with a typical z ^ 2.8 the number of detectable > lOOGeV photons 
will be negligible. Our results are thus largely consistent with the null detection of the MAGIC telescope l| Albert et al. 2007)) . 
Whipple tHoran et al. 2007)) observed the > 400GeV afterglow emission of a few GRBs with 2^1, and in particular GRB 
030329, a nearby long burst. However, the earliest observation was carried out ~ 64.55 hours after the trigger of GRB 030329 
when the expected very high energy SSC afterglow emission is quite low (see also Xue et al. 2008). As the optical depth for 
IR absorption increases strongly with energy and at low redshifts linearly with z, we expect a detection of 100 GeV or lower 
energy photons from nearby strong bursts, provided the observations begin very early and last for several hours. Such nearby 
bursts are, of course, very rare but they do exist and high energy observatories should focus on them. 

7 SUMMARY AND DISCUSSION 

Very high energy inverse Compton emission is an integral part of the current afterglow model l|Meszaros fc Rees 19941 
IDermer et al. 2000|[Sari fc Esin 200T]|Zhang fc Meszaros 2001b|IGou fc MeszaroTDOT. GaUi fc Piro 2007l|Yu, Liu fc Dai 2007] 
IBeloborodov 20051 IFan et al. 2005b] [Wang et al."2006| IFan fc Piran 2006b |l . We have calculated the high energy emission m 
different models of GRB afterglows, including the SSC component of the forward shock, the SSC component of the electrons 
producing X-ray flares and the EIC component of flare photons upscattered by relativistic electrons in the forward shock. 
Our predicted high energy light curves are summarized schematically in the lower panel of Fig. 1151 



* http://wwwmagic.mppmu.mpg.de/ 

^ http:/ /veritas. sao.arizona.edu/old/VERITAS_ whipplc.html 
^ http://www.mpi-hd.mpg.do /hfm/HESS/HESS.html 
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High energy SSC emission in the energy range 20 MeV - 300 GeV from bright bursts should lead to a detectable signal 
of several (« 10) photons by the LAT onboard GLAST. Higher energy telescopes such as MAGIC, Whipple, H.E.S.S and 
Kangooroo working in the energy range > lOOGeV, could detect strong signals (few hundred photons) from nearby bursts at 
around the lower energy limit of these detectors (~ lOOGeV). Signals from more distant bursts will be absorbed by the IR 
background. The flux from a z > 1 burst will usually be too low to be detected. 

Strong GeV SSC emission simultaneous with keV flare photons is possible if the emitting region is not highly magnetized 
0. The EIC component of the flare, on the other hand, will be extended and will last up to ten times as long as the X-ray 
flare (see Fig. I11|FI . This is because, in the EIC process, the duration of the high energy emission is affected by the spherical 
curvature of the blast wave and is mainly extended by the highly anisotropic radiation of the up-scattered photons (see 
Fig |12|l . Unfortunately, a significant detection is likely only if the SSC emission of the forward shock is very weak. A high 
energy detection could be used to probe the spectrum of the low energy flare and in particular the possible existence of a far 
UV component. These signatures of the high energy flare are independent of the density profile of the external medium. 

A detection of a high energy component, m principle, will enable us to test current models of GRBs and their afterglows. 
A detailed comparison of the high energy and low energy light curves, in particular during the shallow decline phase, might 
even enable us to discriminate between different modifications of the standard afterglow model. However, given the small 
number of expected high energy photons, it is unlikely that we can achieve this goal with GLAST. 

It should be noted that the two modifications to the standard model that we considered in this paper, viz., extended 
energy injection and time-evolving e^, both predict achromatic behavior such that there should be a shallow light curve in the 
optical band simultaneously with the shallow X-ray light curve. However, as noted in Fig. 1, this is not always seen. Thus, it 
is possible that none of our models gives a correct description of afterglow physics. Model that suggest fiattening of the light 
curve due to emergence of an X-ray SSC component (see Fig. [3]) that would arise in a very dense medium may explain the 
X-ray flattening but this would involve significant variability in the observed X-ray spectrum. 

One intriguing possibility is that the standard afterglow model, without energy injection or varying e^, is indeed the 
correct model, but the X-ray emission is suppressed during phase-II because of some radiation physics that we have not yet 
understood. This would explain why the optical light curve shows no shallow phase-II segment or a break from phase-II to 
phase-Ill. What kind of high energy emission do we then expect? In the absence of a real model, we cannot say anything 
definite. It is possible that the high energy light curve would follow the predictions of the standard model. Perhaps the high 
energy emission may even be enhanced because the missing X-ray emission is radiated in this band. These are pure speculations, 
but they are worth keeping in mind. It is very important to carry out high energy observations of GRB afterglows, independent 
of model expectations, because the signal may in the end turn out to be stronger than anything predicted. 
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